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The confinement mechanism in the nonperturbative QCD is studied in terms of topological excitation as QCD- 
monopoles and instantons. In the 't Hooft abelian gauge, QCD is reduced into an abelian gauge theory with 
monopoles, and the QCD vacuum can be regarded as the dual superconductor with monopole condensation, 
which leads to the dual Higgs mechanism. The monopole-current theory extracted from QCD is found to have 
essential features of confinement. We find also close relation between monopoles and instantons using the lattice 
QCD. In this framework, the lowest ++ glueball (1.5 ~ 1.7GeV) can be identified as the QCD-monopole or the 
dual Higgs particle. 



1. Dual Higgs Theory for NP-QCD 

Quantum chromodynamics (QCD) is estab- 
lished as the strong-interaction sector in the Stan- 
dard Model, and the perturbative QCD provides 
the powerful and systematic method in analyzing 
high-energy experimental data. However, QCD is 
a 'black box' in the infrared region still now ow- 
ing to the strong-coupling nature, although there 
appear rich phenomena as color confinement, dy- 
namical chiral- symmetry breaking and topologi- 
cal excitation in the nonperturbative QCD (NP- 
QCD). In particular, confinement is the most out- 
standing feature in NP-QCD, and to understand 
the confinement mechanism is a central issue in 
hadron physics. 

In 1974, Nambu [1] presented an interesting 
idea that quark confinement and string picture 
for hadrons can be interpreted as the squeezing of 
the color-electric flux by the dual Meissner effect, 
which is similar to formation of the Abrikosov 
vortex in the type-II superconductor. This dual 
superconductor picture for the NP-QCD vacuum 
is based on the duality in the Maxwell equa- 
tion, and needs condensation of color-magnetic 
monopoles, which is the dual version of electric- 
charge (Cooper-pair) condensation in the super- 
conductivity. 

In 1981, 't Hooft [2] pointed out that color- 
magnetic monopoles appear in QCD as topolog- 
ical excitation in the abelian gauge [2,3], which 
diagonalizes a gauge-dependent variable X(s). 



Here, SU(jV c ) gauge degrees of freedom is par- 
tially fixed except for the maximal torus sub- 
group U(l) Ar<:_1 and the Weyl group. In the 
abelian gauge, QCD is reduced into a U(l) Ar<:_1 - 
gauge theory, and monopoles with unit mag- 
netic charge appear at hedgehog-like configura- 
tions according to the nontrivial homotopy group, 
H 2 {SU(iV c )/U(l)^-i} = Z £-i [4-6]. 

In 90's, the Monte Carlo simulation based on 
the lattice QCD becomes a powerful tool for the 
analysis of the confinement mechanism using the 
maximally abelian (MA) gauge [7-15], which is a 
special abelian gauge minimizing the off-diagonal 
components of the gluon field. Recent lattice 
studies with MA gauge have indicated monopole 
condensation in the NP-QCD vacuum [7-9] and 
the relevant role of abelian degrees of freedom, 
abelian dominance [9-12], for NP-QCD. In the 
lattice QCD in MA gauge, monopole dominance 
for NP-QCD is also observed as the essential role 
of QCD-monopoles for the linear quark potential 
[10], chiral symmetry breaking [11,12] and instan- 
tons [6,13,14]. 

In this paper, we study QCD-monopoles in 
the NP-QCD vacuum using the SU(2) lattice 
QCD in MA gauge. Next, we study the role of 
QCD-monopoles to quark confinement using the 
monopole-current theory extracted from QCD 
[16]. Finally, we study the correlation between 
instantons and QCD-monopoles in terms of re- 
maining nonabelian nature in MA gauge. 
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2. Lattice QCD in MA Gauge 

The maximally abelian (MA) gauge is the best 
abelian gauge for the dual superconductor pic- 
ture for NP-QCD. In this section, we consider 
the mathematical structure of MA gauge. In the 
SU(2) lattice formalism, MA gauge is defined so 
as to maximize 

R = ^{Ms)rsU-\s)r 3 } 

8,fX 

= 2^^? + Ul{sf -Ul{sf -Ul(sf} 

S.fl 

(1) 

by the gauge transformation. Here, = 
exp{iae^l Ai (s)} = U°(s) + iT a U^(s) denotes the 
link-variable on the lattice with spacing a. In 
MA gauge in the lattice formalism, 

4 

A[LV( S )]^]>> ±A1 ( S )T3[/ ± >) 
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(2) 

is diagonalized. Here, we use the convenient no- 
tation as U-n(s) = U~ x {s — jj). 

In MA gauge, there remain U(l)3-gauge sym- 
metry and global Weyl symmetry [13], because 
R is invariant under the gauge transformation 
U^{s) — ► v(s)U l _ l (s)v~ 1 (s + /}) with v(s) = 
e ir 3 ^{s) £ tj^^ and the Wcyl transformation 
U^s) -» WU^W- 1 with W G Wcyl 2 ~ Z 2 
being s-independent. W is expressed as 

W = exp{i7r(— cos</> + — sm(f>)} 

( e~ i ^\ 
= i(TiCos0 + T 2 sin0) = i I 40 Q I (3) 

and interchanges SU(2)-quark color, |+) = (J) 
and |-) = (?). In the SU(A C ) case, this Weyl 
symmetry Weybv c corresponds to the permuta- 
tion group Pm c = Zn c (n c -i)/2, whose element in- 
terchanges SU(-/V C ) -quark color [13,17]. 

Nonabelian gauge symmetry G = SU(A c ) loca i 
is reduced into H = U(l)^ 1 x WeylfJ° bal in MA 



gauge. Then, the independent set of the gauge 
function Oma(s) which realizes MA gauge fixing 
corresponds to the coset space G/H: Oma(s) G 
G/H. The representative element for the link- 
variable in MA gauge is expressed as Ujf A (s) = 
^MA(s)f/ M (s)Oy A (s + fi), and also forms G/H. 
Thus, MA gauge fixing obeys the nonlinear rep- 
resentation on coset space G/H. 

The MA gauge function Oma(s) G G/H is 
transformed nonlinearly by g(s) G G = SU(iV c ) 
as Oma(s) -> ^ma(s) = h[g](s)Q M A(s)g~ 1 (s), 
where h[g](s) G H appears so as to satisfy 
^ma( s ) e G/H. Actually, the successive gauge 
transformation, fl MA after g, is equivalent to 
h[g]QMA, and maximizes R. 

According to the nonlinear transformation in 
JIma G G/H, any operator Oma defined in MA 
gauge transforms nonlinearly as Oma — ► O ma by 
the SU(_/V c )-gauge transformation g G G. 
(Proof) For simplicity, original is assumed to 
obey the adjoint transformation by the SU(7V C )- 
gauge transformation. Then, one finds Oma = 
QmaOQ^ a . By g G G, Oma is transformed as 
Oma - O ma = n 9 MA 63n 9 MA ^ = hlg^MAg- 1 ■ 

gOg- 1 ■ gnuAhig}- 1 = h [g] Oma ^ [5] — ^ = Oma 
with h[g] G H. This proof can be generalized to 
any operator 0. 

If Oma is -H-invariant, one gets M A = ^ma 
for any h[g] G H, so that Oma is invariant under 
arbitrary gauge transformation by g G G. Thus, 
one finds a useful criterion on the SU(A^ c )-gauge 
invariance of the operator in MA gauge [13]. 
"If an operator Oma defined in MA gauge is H- 
invariant, Oma is proved to be also invariant un- 
der the whole gauge transformation of G." 

3. Abelian/Monopole Projection 

The SU(2) link- variable U^(s) can be factor- 
ized as f/ M (s) = M ll (s)u^(s), where u^(s) = 
exp{zr 3 6^(s)} G U(l)3 is abelian link-variable 

and M^s) = e^iW+^M G SU(2)/U(1) 3 , 



/ cos(9 M (s) - e -^( s ) sm6,,{s) 
— y e txn(s) s'mB^is) cos# M (s) 
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with -it < ^(s),Xm( s ) < tt and < 6» M (s) < f . 
Here, cos# M (s) in the abelian gauge is a gauge- 
invariant quantity which measures the 'U(l)- 
ratio' of the link- variable U^(s). For instance, 
(cos^(s)) = 1 means perfectly abelian system. 

In MA gauge, the off-diagonal component of 
Ujf A (s) is strongly suppressed as M^{s) ~ 1 or 
Uf A {s) ~ u^s). Actually, the SU(2) lattice 
QCD in MA gauge shows high 'U(l)-ratio' as 
(cos0 m (s))ma > 0.9 even in the strong-coupling 
region. Therefore, QCD in MA gauge becomes 
similar to the abelian gauge theory. [] 

For any operator 0[Z7„(s)], abelian projection 
is realized as {0\U^(s)]) —> (0[m m (s)])ma- In 
case of (0[C/ M (s)]) ~ (0[m m (s)])ma, the abelian 
degrees of freedom is relevant for 0\U^{s)] in MA 
gauge, which is called as abelian dominance for 
O. For instance, the SU(2) lattice QCD shows 
abelian dominance [10] for the string tension as 
(^K(s)]}MA-0.92.(a[C/ M ( S )]). 

In U(l)3 link- variable u^(s) = exp{ir 3 6'^(s)}, 
d^(s) e (— 7r,7r] is the abelian gauge field on the 
lattice, and abelian field strength is defined as 

0™( s ) ee mod 27r (3 A 6 3 )^(s) G (-7T, tt], (5) 

which is U(l)3-gauge invariant. Generally, 0„(s) 
satisfies (9 A # 3 ) M „(s) = 0^(s) + 2im MV (s). Here, 
n lJLV {s) G Z corresponds to the Dirac string and 
varies by singular U(l)3-gauge transformation. 
There appear magnetic-monopole currents 

K(s) ee ^6l%s) = -d„h MV (a) G Z (6) 

and the electric current j„(s) ee -^d^O^^s). 

We show in Fig.l the monopole current fcu(s) 
in the lattice QCD in MA gauge. In the decon- 
finemcnt phase, monopole currents only appear 
as short-range fluctuation. In the confinement 
phase, monopole currents cover the whole lattice 
and form a global structure, which is an evidence 
of monopole condensation [7,8]. 

Nonperturbative phenomena like confinement 
are brought by large fluctuation of gauge fields 
in the strong-coupling region. In MA gauge, 
such large fluctuation is concentrated into the 

1 In the UV region, such a criterion would be meaningless 
because t/«(s) ^las/3-iooina suitable gauge. 



U(l)3 sector, u fl (s). In particular, monopoles 
appear at the ends of the Dirac strings, and 
accompany large fluctuation of u^(s) or 3 (s). 
Hence, monopole density pu ee y J2 S M I^m( s )I 
is expected to measure the magnitude of gauge- 
field fluctuation. Here, |fc p (s)| and pm in MA 
gauge are SU(A^ C ) -gauge invariant, since |fc M (s)| 
is U(l)3-gauge invariant and \Veyl2 -invariant. 

The abelian gauge field 6n(s) can be decom- 
posed into the monopole part 0ff°(s) and the pho- 
ton part 9 ph (s), 

9™°{S) EE27T^<s|d-V)^( S '), 

=J2A*\d- 2 WWlUs'y (7) 

In the Landau gauge <9 At # 3 t (s) — 0, one finds 
03(a) = 9f l Io (s) + 9 Ph (s). U(l) 3 link-variables 
are defined as uff°> Ph (s) = cxp{iT 3 9ff°' Ph (s)}. 

^From 9ff°(s) and 9 Ph (s), one can derive the 
field strength and the currents in the monopole 
and photon sectors using Eqs.(4) and (5) [11,13]. 
The monopole sector holds the monopole current 
only : fc^ /o (s) ~ fc M (s) and jff°(s) ~ 0, while 
the photon sector holds the electric current only : 
j Ph (s)^Ms) andfc™( S )~0. 

Monopole projection is realized as (0[C/ M (s)]) 
— > (0[u^ Io (s)])ma, and monopole dominance as 
(Op^s)}) ~ (d[uff°(s)}) M A is observed for NP- 
QCD. For instance, monopole dominance for the 
string tension [10] is observed in the lattice QCD 
as (a[uff°(s)}) M A 0.88 • (<r[^(s)]). 

4. Monopole Dynamics for Confinement 

In MA gauge, QCD-monopoles seem essential 
degrees of freedom for NP-QCD. In this section, 
we investigate monopole dynamics and confine- 
ment properties using the monopole-current ac- 
tion [3] extracted from the lattice QCD [16], 

Z= J2 exp{-aJ2kl(s)}6(d^(s)), (8) 

which is defined on lattices with large spacing a. 
In the dual Higgs phase, nonlocal interactions be- 
tween the monopole current would vanish effec- 
tively due to the screening effect [3]. 

Here, the monopole current with length L is re- 
garded as L-step self-avoiding random walk with 
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2c? — 1=7 possible direction in each step. Hence, 
the partition function (7) is approximated as Z = 

Y.LP( L ) e ~ aL - Ei e ~ (Q ~ ln7)L : where the con- 
figuration number of monopole loop with length 
L is estimated as p(L) ~ 7 L . In this system, 
the Kosterlitz-Thouless-type transition occurs at 
a c — In 7 similarly in vortex dynamics in the 2- 
dimensional superconductor. 

We perform direct simulations of partition 
function (7) on lattices. Fig. 2 shows monopole 
density pu = yJ2 s ^\^^( s )\ ano - tnc cluster- 
ing parameter n = Y^i L 2 / Li) 2 as functions 
of self-energy a. As a increases, pu decreases 
monotonously, and declustering of monopole cur- 
rent is observed around a c = 1.8 ln7. 

Now, we study confinement in the monopolc- 
current system using the dual field formalism 
[4,5,17]. We introduce the dual gauge field £? M 
satisfying F^ v = (d A B)^ v . In the dual Landau 
gauge d^B^ = 0, one finds d 2 B^ = k^. Hence, 
starting from the monopole current configuration 
k fJ/ (x), the dual gauge field is derived as 

B,(x) = d-%( X ) =~f d 4 Vj^, 0) 

which leads and the Wilson loop. The Wil- 
son loop (W) shown in Fig. 3 obeys the area law. 
We show in Fig. 4 the string tension aa 2 as the 
function of a. Similar behavior is found between 
Pm and era 2 , which suggests the relevant role of 
monopolcs for confinement. 

Thus, the monopole theory (7) seems to have 
essence of NP-QCD in the infrared region. In real 
QCD, however, the QCD-monopole would have 
its intrinsic size R ~ 0.3fm [3], for it is a collective 
mode composed by gluons. Hence, monopole size 
effects should appear in the UV region, a < R, 
and the monopole action (7) is modified to be 
nonlocal. In fact, the monopole size R may pro- 
vide a critical scale for NP-QCD in term of the 
dual Higgs theory. 

5. Instantons and QCD-monopoles 

The instanton is another relevant topological 
object in QCD according to n 3 (SU(A c )) =Z OQ . 
Recent studies reveal close relation between in- 
stantons and QCD-monopoles [6,8,13-15,18,19]. 



In Fig. 5, we show the lattice QCD result for the 
linear correlation between the total monopole- 
loop length L and Iq = J d i x\tv{G llu G l _ lu )\, 
which corresponds to the total number of in- 
stantons and anti-instantons. The lattice QCD 
shows also monopole dominance for instan- 
tons [13,14]: (Iq[U^) ~ (/g[M^°]) MA and 
(7q[M ai u^ >/i ])ma — after several cooling. 
Hence, instantons can be regarded as 'seeds' of 
QCD-monopoles [6,8,10,13,14,18,19]. 

For these correlation, off-diagonal elements 
would be essential. On the SU(2) lattice in MA 
gauge, we find relatively large off-diagonal ele- 
ments remaining around monopolcs. Hence, in- 
stantons, which need full SU(2) components, ap- 
pear near monopole world-lines in MA gauge. 

The authors would like to thank Professors 
Y. Nambu and R. Brout for their useful comments 
and discussions. 
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FIGURE CAPTION 

Fig.l. Monopole current in MA gauge ex- 
tracted from SU(2) lattice QCD with 16 3 x 4. (a) 
confinement phase (/3 = 2.2), (b) deconfinement 
phase (p = 2.4). 

Fig. 2. (a) Monopole density pu v.s. self- 
energy a. (b) The clustering parameter r\ v.s. 
a. For a < a c , almost monopole world-lines com- 
bine into one large cluster (n ~ 1). For a > a c , 
only small monopole loops appear (n ~ 0). 

Fig.3. The Wilson loop (W(I, J)) v.s. I x J in 
the monopole theory with a=1.7, 1.8, 1.9. 

Fig. 4. The string tension era 2 in the monopole- 
current theory. The dotted line denotes the 
Creutz ratio in the lattice QCD with ,3=1. 25a. 

Fig. 5. Correlation between Iq and the to- 
tal monopole-loop length L in the SU(2) lattice 
QCD. We plot the data at 3 cooling sweep on the 
16 3 x 4 lattice with various (3. 



